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ABSTRACT 


i  i 


The    o^,   O    ,    *o     tracker  system  of  optimal  smoothing  was  investi- 
gated in  a  track-while-scan  radar  by  using  three  dimensional  simulation  of 
various  interceptor-target  profiles.    Random  noise  was  put  on  all  parameters 
which  were  not  positively  defined.    This  included  the  target's  motion. 

The  computer  program  used  is  designed  for  carrying  on  the  analysis 
with  other  systems  so  that  their  performance  can  be  compared  with  the 
o*:  ,   xp  /    vi     tracker.    A  number  of  representative  profiles  are  listed 
which  should  provide  a  suitable  basis  for  evaluation  of  any  chosen 
alternate  system. 

The  parameters  used  in  the  analysis  are  similar  to  those  which 
will  be  used  in  the  Hughes  Aircraft  Co.  developed  AWG-9,  Phoenix, 
Aircraft  Missile  Control  System  designed  for  the  F-111B  aircraft. 
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1 .         Introduction 

The      ^    ,     \p     sampled  data  system  has  been  the  most  often 
proposed  means  of  performing  the  task  of  optimum  smoothing  in  a 
track-while-scan  radar.    This  method  is  particularly  attractive  since 
it  is  very  simple.    However,  in  the  light  of  the  high  performance  intercept 
problems  which  may  be  encountered  today,  especially  by  a  high  perfor- 
mance aircraft,  it  was  felt  advisable  to  investigate  the  capabilities  of 
such  a  smoothing  system  in  the  environment  of  air-to-air  intercepts. 

In  this  analysis  the  third  order  system  or    ^<    ,      \p       ,       0 
tracker  was  considered  by  using  a  three  dimensional  simulation  of  various 
interceptor-target  profiles  .    Random  noise  was  placed  on  all  the  para- 
meters which  were  not  positively  defined.    This  included  the  motion  of 
the  target. 

The  computer  program  used  was  designed  for  carrying  on  the  analysis 
with  other  systems  so  that  their  performance  could  be  compared  with  that 
of  the     Ok  ,    Q>      f  ^      tracker. 

The  parameters  used  in  the  analysis  were  similar  to  those  which  will 
be  used  in  the  Hughes  Aircraft  Company  developed  AWG-9,  Phoenix, 
Aircraft  Missile  Control  System.    This  system  is  currently  being  developed 
for  use  in  the  F-111B  aircraft. 


1  of  100 


*    .    .  B  D  39  0  0  0 

0  0  0  0 

O  0  0   0  0  0 

O  0  0  0 

Dili  •>  i»  .i  ,•  a  a  »> 


2.         Basic  System  Description 

The  system  considered  in  the  analysis  was  a  pulse  doppler  radar 
operating  in  a  track-while-scan  mode    TWS .    That  is,  the  radar  was 
able  to  obtain  simultaneously  at  discrete  sampling  intervals,  called 
frame  times,  the  target's  range,  range  rate,  and  bearing,  in  azimuth  and 
elevation,  relative  to  an  interceptor's  antenna.    At  the  same  time  the 
system  was  also  maintaining  a  track  file  on  a  number  of  other  targets 
at  locations  spaced  over  its  scan  volume. 

The  heart  of  a  system  such  as  this  is  the  so-called  "track-while- 
scan  loop. "    Fig.   1  is  a  basic  block  diagram  of  this  loop. 

Prior  to  entering  the  TWS  loop  the  raw  A  to  D  converted  parameters 
are  taken  from  the  output  of  the  radar,  and  the  interceptor's  inertial 
reference  platform  and  air  data  computer.    This  information  is  pre- 
processed  to  remove  anyredundanciescaused  by  picking  up  the  same 
target  on  two  or  more  antenna  sweeps,  called  bars.    Also  the  north, 
east,  down  direction  cosines  are  computed.    These  are  the  direction 
cosines  for  a  north,  east,  down  axis  coordinate  system  which  is  fixed 
to  the  interceptor.    The  axes  are  defined  by  having  the  north  axis  pointing 
toward  true  north.    The  east  axis  results  from  the  cross  product  of  the 
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down  and  north  axis.    The  down  axis  points  towards  the  earth's  local 
vertical,  and  is  positive  downward.    Thus  an  interceptor  flying  north 
would  have  the  north  axis  along  the  line  of  sight,  the  east  axis  out  the 
right  wing  and  the  down  axis  vertically  downward. 

The  data  now  enters  the  "track-while-*can  loop."    It  is  first 
correlated  by  having  the  individual  observations  checked  against  previous 
tracks  on  file  in  memory.    This  may  result  in  several  observations  being 
linked  with  one  track.    These  are  then  associated  by  having  the  observations 
paired  by  a  logic  routine  with  the  track  which  it  most  closely  matches.    The 
next  stage  initiates  new  tracks  or  deletes  old  tracks,  which  no  longer 
associate  after  a  quality  counter,  which  rates  each  track,  degenerates  to 
a  point  such  that  the  track  no  longer  warrants  being  maintained  on  file. 

The  next  section  and  the  one  which  is  investigated  in  this  paper  is 
that  of  smoothing  and  prediction.    The  ability  of  the  system  to  perform  all 
previous  and  following  tasks  rests  on  how  well  the  track  is  smoothed,  i.e.  , 
the  effect  of  noise  removed,  and  how  well  it  predicts  the  target's  future 
motion.    This  section  must  be  designed  to  respond  quickly  to  target 
maneuvers  while  filtering  out  the  false  excursions  of  the  target  caused 
by  noise  so  that  future  predictions  will  be  as  accurate  as  possible. 
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The  gating  section  is  next  in  the  loop.    It  is  this  section  which 
determines  the  possible  limits  of  target  motion  between  contacts.    It  allows 
for  target  maneuvers  between  scans  plus  approximately  3cr noise.    It  further 
provides  for  expansion  in  the  event  of  a  missed  contact.    A  maximum  gate 
size  is  imposed  which  corresponds  to  a  volume  which  would  contain  any 
possible  target  maneuver  that  could  occur  up  to  the  time  the  track  would 
normally  be  eliminated  by  the  quality  counter  due  to  missed  contacts. 

The  individual  track  parameters  are  stored  in  the  track  file  section 
to  keep  it  updated .    These  along  with  the  track  locations  are  used  to 
determine  target  priorities .    The  track  parameters  are  also  used  to  determine 
which  observations  correlate  with  the  tracks  during  the  correlation  phase 
at  the  beginning  of  the  loop.    The  above  computations  are  performed  on  all 
tracks  at  the  end  of  each  frame  time . 

Additional  functions  which  result  from  the  TWS  loop  computations 
are  the  antenna  scan  center  calculations  to  keep  the  pattern  centroid 
located  in  the  center  of  the  incoming  raids ,  the  steering  information  to  the 
pilot  to  allow  him  to  maneuver  as  necessary  to  keep  the  radar  within  its 
gimbal  limits,  prelaunch  information  to  the  missile  so  that  its  seeker  head 
is  oriented  in  the  proper  position  at  the  beginning  of  its  semi-active  phase, 
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and  midcourse  maneuver  corrections  to  the  missile  during  flight. 

It  is  readily  apparent  that  the  system  must  be  optimized  to  the  highest 
possible  degree  in  light  of  the  high  relative  speeds  encountered  and  the 
number  of  targets  to  be  tracked.    At  the  same  time  it  is  extremely  important 
that  computations  remain  relatively  simple  to  maintain  reliability  and  to 
keep  computational  time  and  storage  space  to  a  minimum. 

The  ability  of  the  system  to  meet  these  criteria  lies  primarily  in  the 
anoothing,  prediction  and  gating  sections.  It  is  the  former  two  which  are 
considered  in  this  paper. 
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3.         System  Characteristics 

The  values  of  the  system  performance  parameters  were  derived 
primarily  from  the  Hughes  SYSTEM  EVALUATION  AND  ANALYSIS  REPORT, 
SEAR,  /l^/.    Where  possible  these  parameters  have  been  crossed  checked 
with  other  company  internal  documents ,  or  by  conversations  directly  with 
company  engineers.    The  following  parameters  and  assumptions  are  those 
which  were  used  in  the  simulation. 

It  was  assumed  that  the  interceptor  would  have  an  on  station  Mach 
number  of  0.7  at  above  36,000  feet.    The  target's  speed  was  believed  to 
be  between  Mach  0.4  and  3.0  and  its  altitude  could  be  between  sea  level 
and  80,000  feet. 

The  general  characteristics  of  the  radar  were  those  of  a  high  PRF 
pulse  doppler  using  linear  frequency  modulation. 

The  detection  range  was  165  n.m.    This  was  defined  as  the  range 
where  S/N  was  one  using  the  classic  radar  range  equation.    It  was  based 
on  a  five  square  meter  target.    The  radar's  probability  of  detection  curve 
is  shown  in  Fig .  2  . 

The  effect  of  clutter  on  detection  probabilities  was  not  considered 
since  a  large  amount  of  clutter  rejection  is  an  inherent  advantage  of  a  high 
PRF  pulse  doppler  radar.    This  is  particularly  true  over  water.    In  addition 
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the  system  will  also  have  clutter  rejection  circuitry.    Due  to  the  high  PRF 
certain  relative  velocities  may  cause  periodic  loss  of  signal  due  to 
"eclipising" ,  i.e.  ,  target  echo  returning  during  pulse  transmission. 
Multiple  PRF  rates  will  normally  remove  this  possibility.    Thus  eclipising 
was  not  considered  to  effect  detection  probability. 

In  the  radar,  range  is  measured  by  taking  the  difference  between 
two  filters  in  a  doppler  filter  bank  that  are  "rung"  on  successive  FM  phases 
The  spacing  of  the  filters  corresponds  to  five  n.m.    Since  there  are  30 
filters  the  range  can  be  measured  to  150  n.m.    The  signal  was  considered 
to  originate  uniformly  in  plus  or  minus  one-half  a  filter  and  terminate 
normally  distributed  with  a  standard  deviation  of  one-third  a  filter. 
This  is  shown  in  Fig.   3. 
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Fig.   3 
Probability  Distribution  of  Range  Measurement 
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Hughes  stated  that  a  joint  distribution  describing  these  two  distributions 

yielded  a  normal  distribution  of  zero  mean  and  standard  deviation  of 

1 .  37  n.m.    It  was  felt  that  this  value  was  too  optimistic  and  indicated  that 

the  distributions  were  quite  dependent.,  i.e.  ,  the  covariance  was  large. 

The  value  used  for  the  range  standard  deviation  in  this  analysis  was 

2.5  n.m. 

Range  rate  is  measured  by  frequency  shift  as  determined  by  a  doppler 
filter  bank.    The  filters  are  spaced  150  cps  apart.    Thus  at  X  band  using 
the  equation  given  by  Skolnik  /_2_J: 


V        =     d 

r 


f"     >S 


103 


Where:      i>    =    doppler  frequency  shift. 

/\    =    wavelength  in  cm.  ,  3  cm  at  X  band 


V     =    range  rate  in  knots  . 
r 


The  total  number  of  filters  is  768  giving  the  ability  to  measure  range 
rates  of  up  to  approximately  3250  kts . 

It  was  determined  by  Hughes  that  the  distribution  of  the  range  rate 
error  was  normal  with  zero  mean  and  a  standard  deviation  of  two  filters , 
8.72  kts. 
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The  scan  pattern  of  the  antenna  in  the  TWS  mode  is  either  two  bar 
+  40  deg.  or  four  bar  +  20  deg.    The  beam  separation  is  1.32  deg.  and  the 
scan  rate  is  40  deg/sec.    It  was  determined  from  various  sources  that  the 
average  antenna  azimuth  error  could  be  considered  normally  distributed  with 
zero  mean  and  standard  deviation  of  0.745  deg.    Similarly  the  elevation 
error  had  zero  mean  and  standard  deviation  of  0.86  deg.    These  values 
include  radome ,  potentiometer  output,  and  pick  off  errors.    The  additional 
error  in  elevation  is  due  primarily  to  the  introduction  of  redundancies  caused 
by  picking  up  the  same  target  on  two  separate  scans. 

The  smoothing  and  prediction  system  which  has  been  proposed  by 
Hughes  is  the    c^  t    (X     ,    "^     tracker  type.    The  characteristics  of  this 
system  are  considered  in  the  following  section. 
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4.         Smoothing  System  Discussion 

Since  computational  time  and  storage  space  available  for  smoothing 
and  prediction  is  quite  limited  some  of  the  more  sophisticated  smoothing 
systems  would  not  at  first  appear  to  be  suitable  for  use  in  the  system 
described.    The  requirements  of  the  system  are  reduced  somewhat,  however, 
by  the  fact  that  the  target's  capabilities  are  relatively  predictable,  but  of 
course,  entirely  random.    For  these  reasons  the  most  attractive  choice  is 
an     <=x    ,  9}      tracker  type  of  system.    Almost  any  other  system  that  would 
be  acceptable  would  reduce  to  a  scheme  which  would  use  a  similar  smoothing 
constant  gain  technique  in  any  case. 

The    "^   ,  \p     o  "^Q     tracker  system  is  an  extension  of  the      ^X  ,    P 
tracker  to  allow  for  the  smoothing  of  acceleration .    As  given  by  Benedict 
and  Bordner  /  3/  the  variance  reduction  ratio  is  optimized  by  having  the 
velocity  constant,    \:b     ,  defined  as  : 


^> 


^*r 


1  -•< 

Similarly  it  was  shown  by  Simpson  /4/  that  when  acceleration  was 
considered  the  optimum  value  of  the  acceleration  smoothing  constant,   o"  , 


became: 


y-^-.h^. 
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Thus  if  it  is  possible  to  define  ^><    the  smoothing  equations 
become  unique . 

The  proposed  equations  used  for  the  smoothing  and  prediction 


are: 


A.      Range: 


♦£"*-' 


ta^s    -Rn   *    T"?.. 


~V» 


Where: 


-R.J*- 


o^  and 


IV-  h       V^  o  ^        i 


If  an  observation 


was  received 


l<Ln=-  f^p      and        K^  ^  -   ^9^       if  an  observation  was  missed 

f- Frame    Time    =    2  +0.1  sec. 

cj  =  Number  of  frame  times  since  last  value  of  range  rate  was 
received. 

The  subscripts  are: 
P>  -  Frame  number. 
P=  Predicted. 
O  =  Observed. 
<,-  Smoothed,  I  % 
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The  value  of    ^    for  range  smoothing  as  proposed  by  Hughes  was 
0.45.    Since  the  standard  deviation  of  range  rate  was  relatively  small 
no  smoothing  of  range  rate  was  necessary  so  that  no    O     ,  or        0 
were  required  in  range  smoothing. 

B.      Direction  Cosines: 


Jv,=  A^-o  ^  (Jw  JVP^ 


Where: 

__7V  ns    -J\.0t.    if  an  observation  was  received. 

/\       w  J\_p      if  an  observation  was  missed. 

_J\_    =      general  direction  cosine  being  smoothed. 
For  the  North  and  East  direction  cosines,  N,  E,  Hughes  proposed 
the  smoothing  constant  values  of: 
<=><  ~     0.43 

^X  -^  0.120 
.^-    0.008. 
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For  the  down  direction  cosine,  D,  the  proposed  values  were: 
^s~-  0.35 
&  =0.075 
^-0.0035. 

A  Z-  transform  analysis  of  the  c^  ,    P>      Tracker  was  performed  by 
Sklansky  /5_/.    Basically  the  system  is  stable  and  for  a  non-maneuvering 
target  there  would  be  no  steady  state  error  for  ^<  between  zero  and  one. 
When   <:=*v  is  equal  to  one  no  smoothing  would  be  performed  and  raw  data 
would  be  the  output.    In  a  high  noise  non-maneuvering  target  condition 
the  smoothing  constant  would  be  small  and  the  system  response  would 
be  slow.    For  increasingly  lower  noise  the  smoothing  constant  would 
approach  one  in  the  limit.    For  a  maneuvering  target  the  smoothing 
response  must  be  increased  so  that  a  larger  value  of  Cs^  is  required. 

The  obvious  disadvantage  of  the  ^^  ,  \p    Tracker  is  that  there  is 
not  a  single  optimum  value  of  the  smoothing  constant  for  all  the  situations 
that  can  be  encountered. 

The  advantage  of  the  system  again  lies  in  its  simplicity  and  the 
small  number  of  computations  and  amount  of  storage  space  required. 
Also,  since  o    and  ^     are  defined  uniquely  when  "^    is  determined  the 
problem  of  fixing  the  constants  is  reduced  somewhat. 
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5 .         Simulation 

The  CDC  1604  computer  was  used  to  implement  the  three  dimensional 
simulation  of  a  single  target-interceptor  problem.    This  program  was  named 
Three  D.    The  primary  aims  of  the  simulation  were  to  provide  realistic, 
flexible  and  accurate  target-interceptor  geometries .    The  pseudo  random 
noise  generator ,  sub  routine  RNDEV ,  was  used  to  provide  gaussian  noise 
on  all  parameters  that  were  not  positively  defined. 

A  list  of  symbols  and  a  description  of  the  use  of  the  programs  used 
are  contained  in  Appendices  A  and  B. 

The  resulting  program  allowed  the  target  to  be  positioned  initially 
at  any  location  in  the  northern  half  of  a  north,  east,  down  coordinate 
system  which  was  ground  stabilized  at  the  interceptor's  initial  location 
at  problem  time  zero.    The  target's  initial  coordinates,  maneuver  para- 
meters, and  maneuver  times  were  then  read  in  by  data  cards  along  with 
the  interceptor's  heading,  altitude,  and  Mach  number.    The  interceptor's 
velocity  was  computed  in  the  program  by  using  the  following  relation: 

Vfps    =    Sonic  Velocity  fps  x  Mach  Number. 
The  variation  of  the  speed  of  sound  with  altitude  is  shown  in  Fig.  4. 
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Speed  of  Sound  Versus  Altitude 


No  provision  was  included  to  maneuver  the  interceptor  since  one  of 
the  basic  assumptions  of  the  system  was  that  the  interceptor  was  continually 
tracking  several  targets  in  widely  spaced  locations  in  the  radar' s  scan 
pattern.    It  was  thus  not  possible  for  the  pilot  to  place  an  individual  target 
in  the  center  of  the  scan.    Also  no  random  noise  was  used  to  vary  the  inter- 
ceptor's position  since  the  same  relative  effect  was  gained  by  randomly 
varying  the  target's  velocity  components. 


17  of  .100 


»  >   *  I   *  11  >   .   1  >   *  >   >   *  *   * 

<  >  >     >  X     1 

x     .     > 


The  problem  time  was  taken  to  be  300  seconds.   This  allowed  a  target 
at  Mach  1.5  above  36,000  feet  to  travel  slightly  over  70  n.m.  in  the  course 
of  the  problem.     Since  the  probability  of  acquiring  and  tracking  a  target 
beyond  100  n.m.  was  very  low,   100  miles  could  be  considered  the  range 
limit  in  setting  up  the  problem  profile. 

The  target's  velocity  was  also  determined  by  the  previously  mentioned 
relations.    No  provision  was  made  to  accelerate  the  target,  but  the  target 
would  decelerate  at  the  start  of  its  dive.    Its  final  speed  was  a  function 
of  altitude  below  10,000  feet.    The  minimum  speed  was  Mach  1.2,   1340  fps. 
at  sea  level.    The  relations  used  are  shown  in  Fig.  5.     The  target  was 
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programmed  to  make  a  180  degree  turn  at  a  programmed  time  and  horizontal 
turn  g  load.    It  would  also  make  a  60  degree  dive  and  pull  out  at  a  pro- 
grammed time,  final  altitude,  positive  and  negative  g  factor. 

For  convenience,  distances  were  read  in,  in  nautical  miles,  velocities 
in  Mach  number,  and  headings  in  degrees.    The  program  then  converted  the 
values  to  compatible  units  of  ft.  ,  ft. /sec.  ,  and  radians. 

The  motion  of  the  target  was  computed  by  piece-wise-linear  integration. 
Random  accelerations  having  a  standard  deviation  of  one-half    g    for  two  second 
periods  were  imposed.    The  target's  motion  was  broken  down  into  velocity, 
VT ,  along  its  path  of  flight,  heading  angle,  AT,  and  pitch  angle,  BT.    Each 
component  was  then  related  to  the  distance  traveled  in  each  of  the  coordinate 

axes  in  the  iteration  time  and  combined  with  the  previous  position  to  give 

i 

the  new  target  location  ,  XT  ,  YT  ,  ZT  . 

Every  two  seconds  of  problem  time  the  target's  coordinates  were  taken, 
compared  with  the  interceptor's  coordinates  XO,  YO ,  ZO  and  the  actual 
antenna  azimuth,  ATA,  elevation,  EPSI ,  range,  R,  and  range  rate,  RDOT 
were  computed.    If  the  target  had  passed  out  of  the  limits  of  the  radar  beam,+ 
65  degs.  in  azimuth  and  elevation  the  problem  stopped.    It  should  be  noted 
that  the  antenna  does  not  have  a  scan  volume  this  large  in  the  track  while 
scan  mode ,  but  it  does  have  the  capability  of  positioning  its  scan  anywhere 
inside  these  limits 
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The  random  noise  generator  was  then  used  to  put  noise  on  the  radar 
parameters,  ATAE ,  EPSIE ,  RE,  and  RDOTE ,  and  the  apparent  target  location 
having  coordinates  XTE ,  YTE,  and  ZTE ,  was  determined.    Since  a  relatively 
small  number  of  samples  were  used  the  mean,  variance  and  independence 
of  the  noise  for  each  of  the  parameters  were  checked  using  an  autocorrelation 
routine.    The  results  of  analysis  are  shown  in  Table  I. 

The  actual  and  noisy  direction  cosines  for  the  three  axes'  were  then 
computed.    These  were,  north  direction    cosine,  DIRN ,  and  DIRNE ,  east 
direction  cosine,  DIRE  and  DIREE,  and  down  direction  cosine,  DIRD,  and 
DIRDE. 

Knowing  the  probability  of  detection  of  the  radar,  RNDEV  was  used 
to  determine  a  zero-one  probability  of  detection  for  each  frame  time.    This 
was  done  by  cross  plotting  the  absolute  value  of  the  deviation  required, 
versus  range,  corresponding  to  the  probability  of  detection  curve,  Fig.  2. 
This  was  then  made  piece  wise  linear  and  implemented  by  using  IF  statements 
Thus  for  each  actual  range  a  maximum  value  of  deviation  for  detection  could 
be  determined,  if  the  absolute  value  of  DEV  obtained  from  RNDEV  exceeded 
this  value  no  return  for  that  frame  time  was  assumed.    Fig.  6  illustrates  the 
method  used. 
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TABLE  I 

AUTOCORRELATION  ANALYSIS  RESULTS  OF  PROGRAM 

THREE  D 
(151  SAMPLES) 


Mean 

Variance 

Largest  Value  of 
Autocorrelation 

RN  1 

Target 

Acceleration 

0.0112 

.772 

0.1787 

Noise 

RN  2 

E5 

Target 

Heading 

Acceleration 

-0.001 

1,0636 

0.1706 

Noise 

RN  3 

Target 

Pitch 

Acceleration 

-0.1032 

0.9585 

0.2358 

Noise 

ATA 

Antenna 

Azimuth 

0.0785 

0.9612 

0.1717 

Noise 

EPSI 

Antenna 

Elevation 

0.6308 

1.0433 

0.1576 

Noise 

R 

Radar 

Range 

0.1386 

0.9938 

0.1721 

Noise 

RDOT 

Radar 

Range 

Rate 

0.1114 

0.9780 

-0.1528 

Noise        

•               '                   1 

PD 

t           0 

! 

•'!         i 

Radar 

...           J  9 

• 

»                   .» 

Probability 

of  Detection 

-0.0564 

0.9914 

-0.1468 

Noise 
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In  the  actual  system  this  would  be  the  status  of  the  individual  track 
in  the  TWS  loop  after  being  correlated  and  associated. 

A  quality  counter,  NQ,  was  devised  for  the  program  which  incremented 
in  steps  of  two  for  each  received  return  to  a  maximum  of  ten.    For  each 
missed  return  it  degenerated  by  one.    When  the  quality  counter  reached 
zero  the  track  was  dropped  and  no  smoothing  was  performed  until  the  next 
received  return. 

The  particulars  of  the  smoothing  and  prediction  section  were  as 
described  earlier. 

The  primary  quality  indicator  for  the  system  was  the  three  dimensional 
rms  smoothed  difference  from  the  actual  position.    Since  high  relative 
velocities  between  the  missile  and  the  target  are  considered  very  likely 
in  the  active  phase  of  the  missile  flight  it  is  necessary  that  the  missile  be 
positioned  accurately  at  the  beginning  of  this  period  so  that  its  maneuvers 
may  be  kept  to  a  minimum.    In  order  to  assist  in  achieving  this  a  desired 
rms  smoothing  error  limitation  of  approximately  one-half  mile  inside  50  miles 
was  imposed.    On  the  output  curves  points  were  plotted  indicating  these 
limits.    For  range  constant  lines  of  points  at  plus  or  minus  one-half  mile 
were  plotted.    For  azimuth  and  elevation  points  were  plotted  as  a  function 
of  range  using  the  relation: 
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If  all  the  parameters  were  smoothed  to  the  desired  value  the  three  dimensional 
rms  error,  CIRSMO,  could  be  as  high  as  0.875n.m.    Thus  on  the  curves  of 
CIRSMO  a  line  of  points  having  this  value  were  plotted. 

The  program  print  out  contained  the  run  number,  problem  time,  actual 
range,  rms  three  dimensional  actual  error,  CIRERR,  rms  three  dimensional 
smoothed  error,  CIRSMO,  a  one  or  zero  for  detection  or  no  detection,  ND, 
and  the  status  of  the  quality  counter,  NQ.    When  the  quality  counter  was 
zero  no  contact  was  assumed  and  all  nines  were  printed  out. 

To  observe  the  results  of  the  smoothing  method  the  sub-routine 
CALL  DRAW  was  used  to  plot  the  interceptor  and  target  positions  in  the 
north-east,  N-E,  plane,  looking  down  from  above,  and  the  north-down, 
N-D,  plane,  looking  in  from  the  east  axis.    The  subroutine  was  also  used 
to  plot  the  difference  between  the  actual  and  smoothed  values  of  range 
and  the  direction  cosines,  which  were  plotted  as  angular  differences,  and 
the  values  of  the  three  dimensional  rms  smoothed  error. 
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6.         Results  of  the  Analysis  of  the  Proposed  System 

Table  II  contains  the  important  details  of  the  problem  profiles  used 
in  the  analysis.  The  results  of  run  1A  are  shown  in  Table  III  and  Figures 
7  to  15.    The  results  of  runs  2A  and  3A  are  shown  in  Appendix  D. 

Three  basic  profiles  were  considered.    These  were  identified  by  the 
initial  position  and  heading  of  the  target  with  respect  to  the  interceptor. 
For  run  1  the  target  was  positioned  at  100  n.m.  and  was  closing  the 
interceptor.    For  run  2  the  target  was  positioned  at  10  n.m.  and  was 
moving  away  from  the  interceptor.    In  run  3  the  target  was  positioned  at 
approximately  60  n.m.  at  a  30  degree  angle  from  the  interceptor  and  was    . 
headed  on  a  cross  track  with  respect  to  it.    The  actual  profiles  flown  by 
the  target  and  the  interceptor  could  be  varied  in  a  myriad  of  ways ,  however 
in  reaching  the  conclusions  that  follow  only  those  flight  paths  shown  in 
Table  II  were  used. 

The  results  obtained  for  the  system,  as  proposed  by  Hughes,  indicated 
that  the  smoothed  tracks  were  still  quite  noisy.    The  direction  cosine  smoothing 
results  obtained  by  analyzing  the  figures  of  these  values  compared  to  the 
desired  limits  indicated  that  the  direction  cosines  were  being  smoothed  inside 
the  limits  rather  well.    Some  degradation  was  noted  in  the  cross  track  cases 
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SMOOTHING  RESULTS 

RUN  1A 

TsJlME*  R  NM=ACTUAL  RANGE,  C I RERR  =  3D I  MENS IONAL  RMS  ACTUAL  ERROR, 
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where  the  rates  of  change  of  the  direction  cosines  were  high.    As  was 
expected  the  range  smoothing  results  were  the  prime  cause  of  the  smoothed 
error . 

The  system  was  able  to  smooth  the  tracks  to  within  one  and  a  half 
nautical  miles  inside  50  miles  rather  consistently,  but  it  tended  to  respond 
rather  strongly  to  the  noise.    It  was  noted,  however,  that  this  had  the 
effect  of  causing  no  degradation  in  performance  when  the  target  began 
maneuvers.    Also  all  the  profiles  yielded  very  similar  smoothing  results 
indicating  that  the  smoothing  constant  values  used  were  near  optimum  for 
the  fixed  constant,  non-adaptive,  system. 

In  an  attempt  to  more  closely  define  the  optimum  values  of  the 
smoothing  constants  a  program  which  varied  the  constants  for  each  smoothed 
parameter  from  0.05  to  1.00  in  steps  of  0.05  was  written.    This  program 
is  named  ALPHA  CHECK  and  listed  in  Appendix  C .    The  run  profiles  used 
were  similar  to  those  listed  earlier    but,  the  run  time  was  shortened  to 
100  seconds  so  that  only  the  time  that  the  target  maneuvered  was  considered, 
The  average  mean  square  error  for  the  complete  run  was  computed  for  each 
smoothed  quantity  over  the  range  of  the  values  of  its  smoothing  constant. 
An  autocorrelation  analysis  was  also  performed  to  evaluate  this  noise 
distribution...   Thooror.'alts  are.  shov^n  ?.n  Table  IV,,,    Eor  these  juns  the 
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AUTOCORRELATION  ANALYSIS  RESULTS  OF  PROGRAM 

ALPHA  CHECK 
(51  SAMPLES) 


Mean 

Variance 

Largest  Value  of 
Autocorre  lation 

RN  1 

Target 
Acceleration 

-0.1735 

0„6945 

0.3947 

Noise 

RN  2 

Target 

Heading 

Acceleration 

0.2439 

1.1572 

0.2662 

Noise 

RN  3 

Target 

Pitch 

Acceleration 

-0.0986 

.9622 

0.1774 

Noise 

ATA 

Antenna 
Azimuth 

0.264 

1.2022 

-0.2493 

Noise 

EPSI 

Antenna 
Elevation 

-0.0765 

0.8361 

0.2544 

Noise 

R 

Radar 
Range 

0.1438 

0   7  524 

0.2783 

Noise 

RDOT 

Radar 

Range 
Rate 

-0.0661 

1 o 1695 

0.2805 

Noise 

PD 

Radar 

Probability 
of  Detectiqir  ; 

-;•   -0.;05e'0,"- 

f 
0.7587  ■•  •• 

|>      .0.-3582 

Noise 

;    :       ' :    : 
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the  probability  of  detection  was  one.    This  did  not  effect  the  results  since 
no  smoothing  would  have  been  performed  when  a  contact  was  missed.    The 
output  would  have  been  affected  in  an  indeterminable  way,  however,  by 
these  predicted  values. 

The  results  of  these  runs  are  shown  in  Table  V.    It  was  apparent 
from  these  results  that  there  was  considerable  interplay  between  all  the 
variables.    Some  rationale  could  be  seen  in  this  interaction,  but  it  appeared 
nonlinear,  since  any  change  in  the  target's  profile  varied  all  the  points  of 
the  optimum  values  for  the  smoothing  constants.    It  was  noted  further  that 
the  knuckle  in  the  curves  at  the  optimum  point  was  quite  sharp  indicating 
that  system  performance  could  be  seriously  affected  when  improper  values 
of  the  smoothing  constants  were  used. 

Due  to  the  variations  which  were  found  in  the  optimum  values  of  the 
smoothing  constants  for  the  chosen  profiles  it  was  not  believed  to  be 
necessary  to  vary  the  parameters  more  widely  since  it  was  apparent  from 
these  runs  that  the  system  was  very  sensitive  to  changes  in  the  profiles 
or  the  values  of  the  system  noise. 

The  question  must  then  be  asked  as  to  whether  the  requirements 
established  for  the  analysis  were  too  strict.    In  this  regard  specifically 
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OPTIMUM  SMOOTHING  CONSTANT  VALUES 

FOR 
EACH  TEST  RUN 


Run 

Range 

Direction  Cosines 

North 

East 

Down 

1  A 

0.1 

0.25 

0.25 

0.35 

-LB- 

0.05 

0   25 

0.25 

0.05 

LP 

0.15 

0.35 

0.35 

0.05 

2  A 

0.25 

0   90 

0.85 

0.65 

2  B 

0.05 

0.65 

0.55 

0.10 

2  C 

0.35 

0.95 

0.95 

0.40 

3  A 

0.50 

0.60 

0.65 

0.55 

3  B 

0.45 

0   70 

0.55 

0.05 

3  C 

0.55 

0.75 

0.80 

0.05 
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the  standard  deviation  of  range  noise  and  the  0.5  mile  value  chosen  for 
the  desired  smoothing  error  become  suspect  since  they  were  somewhat 
arbitrarily  chosen.    In  rebuttal  it  is  strongly  felt,  in  the  light  of  the  high 
overall  performance  required  of  the  entire  system,  its  cost,  and  purpose, 
that  the  use  of  this  relatively  uncontrollable  method  of  smoothing  could  be 
disastrous.    Since  the  output  of  the  smoothing  section  has  considerable 
control  over  the  performance  of  the  system,  specifically  the  accurate 
positioning  of  the  missile,  it  is  believed  that  a  more  adaptive  system  is 
desirable .    It  is  further  felt  that  a  programmed  logic  system  which  would 
provide  for  varying  the  values  of  the  constants  would  not  be  warranted 
since  more  sophisticated,  and  proven  adaptive  systems  are  available. 
Possible  choices  for  such  a  method  could  be  either  the  Least  Squares , 
or  Kalman  filter  methods  . 

In  short  it  must  be  concluded  that  the     C:V  ,     0}       ,     Q>       Tracker 
system  does  not  have  the  response  or  versatility  required  for  a  high 
performance  air  to  air  intercept  problem. 
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NOTE! 


ADEV: 


APPENDIX  A 


TABLE  OF  SYMBOLS  USED  IN  COMPUTER  PROGRAMS 


If  a  symbol  is  used  only  in  the  ALPHA  CHECK  program  it  is  so 


noted  in  parenthesis . 


Absolute  value  of  deviation  obtained  from  RNDEV. 


This  is  used  in  the  probability  of  detection  computations 


ADSMO: 
AESMO: 
ANSMO: 


Down,  east,  and  north  direction  cosine  differences 


from  actual,  and  converted  to  radians. 


AD: 
AE: 
AN: 
AR: 


(ALPHA  CHECK)    Values  of  alphas  for  down,  east,  and 
north  direction  cosines,  and  for  range,  used  for  determining 
the  optimum  values . 


AERRM:  Negative  and  positive  desired  limits  on  ADSMO,  AESMO, 

AERRP:  and  ANSMO,  a  function  of  actual  range. 
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ALP  HAD: 
ALPHAE. 
ALPHAN: 
ALPHAR: 


Position  smoothing  constants  for  down,  east  and  north 
direction  cosines,  and  for  range. 


ANOISE: 


(ALPHA  CHECK)    Relative  amount  of  radar  parameter  noise 


used  in  problem. 


AT: 
ATDOT: 

ATA: 
ATAE:  . 
AT  AD: 
ATADE: 


Target  heading  and  turn  rate  used  in  computation  of  target's 


position. 


Antennas  actual  and  noisy  azimuth  angle  measured  positive 


clockwise,  in  radians  and  degrees. 


BET  AD: 
BETAE: 
BET  AN: 


Rate  smoothing  constants  for  down,  east,  and  north 
direction  cosines. 


BT: 
BTDOT: 


Targets  pitch  angle  and  pitch  angle  rate  used  in  computation 


of  targets  position. 
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CIRERR*  Three  dimensional  RMS  actual  and  smoothed  error  difference 

CIRSMO:  from  actual  target  position. 


DEV: 


Random  deviation  obtained  from  random  number  generator 


DT: 


Target  position  computation  step  size. 


DIRN: 
DIRNE: 
DIRNS: 
DIRNP: 


North  direction  cosine  actual,  noisy,  smoothed  and 


predicted 


DIRE: 
DIREE; 
DIRES: 
DIREP: 


East  direction  cosine  actual,  noisy,  smoothed  and 
predicted. 


DIRD: 
DIRDE: 
DIRDS: 
DIRDP: 


Down  direction  cosine  actual,  noisy  smoothed  and 
predicted . 
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DIRND: 


DIRNDD: 


DIRED: 


DIREDD: 


DIRDD: 


DIRDD: 


North  direction  cosine,  rate  and  acceleration 


East  direction  cosine,  rate  and  acceleration 


Down  direction  cosine,  rate  and  acceleration 


DDSMO:  (ALPHA  CHECK)    Down  east  and  north  direction  cosine 

DESMO:  smoothing  differences  from  actual. 

DNSMO: 

DELTAT:  Target  decceleration  time  to  final  velocity  at  -1.75  g. 


DIVE: 


Target  maximum  dive  angle,  radians,  corresponds  to  60  deg, 


EPSID: 
EPSIDE: 
EPSI: 
EPSIE: 


Antenna  actual  and  noisy  elevation  angle  measured  positive 
downward,  in  radians  and  degrees. 
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ERAD: 
ERAE: 
ERAN: 
ERAR: 

F: 


(ALPHA  CHECK)  Average  mean  square  smoothed  error  for 
down,  east,  and  north  direction  cosines,  and  for  range. 
They  are  the  quality  indicators  used  for  obtaining  the 
optimum  values  of  the  smoothing  constants. 

A  counter  which  fixes  the  target's  pull  out  velocity. 


FPOINTS:  (ALPHA  CHECK)    Number  of  radar  returns  received . 

GAMMAD:  Acceleration  smoothing  constants  for  the  down,  east, 

GAMMAE:  and  north  direction  cosines . 
GAM  MAN: 

GDN*:  Target  maneuver  load  factors  for  push  over,  pull  out, 

GDP:*  and  turn.    Note:    GDN  is  a  negative  number. 

GT:* 


IT: 


Hollerith  characters  which  identify  the  output  plots . 


I: 


Number  of  frame  times  since  problem  start, 


K:  Number  of  desired  limits  points  plotted  from  problem  start 

*    Read  into  the  program 
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NALPH:  (ALPHA  CHECK)    Number  of  values  of  smoothing  constants 

used  in  finding  optimum  value. 


ND: 


A  zero  or  one  depending  on  whether  detection  was  missed 


or  not 


NERMAX:  Number  of  points  plotted  indicating  the  desired  limits  on 

the  smoothed  errors . 

NGRAF1:*         Setting  NGRAFS  to  zero  supresses  various  portions  of  the 

NGRAF2:*         graph  output. 

NGRAF3:* 

NOISE:*  Number  of  times  RNDEV  is  stepped  prior  to  problem  start  to 

allow  a  different  sequence  of  random  numbers  to  be  called. 

NPRINTs  Number  which  separates  pages  of  print  out. 


NQ: 


Quality  counter. 


NQC: 


Counter  which  reinitializes  the  initial  conditions  after  the 


quality  counter  goes  to  zero. 
*    Read  into  the  program 
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NUMPTS:  Number  of  curve  points  plotted. 


NUNIF:  Counter  in  random  number  generator. 


OM 


.* 


Intercept  >r  Mach  number  „ 


PD: 


Target  detection  probability  for  each  frame  time,  a  function 


of  the  actual  range 


PI: 


3.14159265. 


POVERD:  Target  pushover  rate,  rad/sec. 

POUTD:  Target  pull  out  rate,  rad/sec . 

PSIO:  Interceptor  heading  in  radians  and  degrees 
PSIOD:* 


PSIT: 


Target  heading. 


PSITD:*  Target  initial  heading,  deg.  and  rad.    Note  a  zero  not  an 

PSITO:  alphebetic  o  on  target°s  initial  positions. 


Q:  Number  of  frame. times  since  last  radar  return  of  range  rate. 

*    Read  into  the  program   , 


>    i » »  *    >  i  >    ill*    i  i  > 


i      i  >  j      ii>)      i    1 1 1  <    1 1 1      i      i 
•i  j  i  i  i 

I.       jiio  ii,i 


51  oFlOQ 


>    3  i  i 





RE: 

RNM: 

RNME: 


Actual  and  noisy  range  to  target,  ft  and  n.m 


RERRM: 


RERRP: 


Negative  and  positive  desired  limits  on  range  „  a 


constant  .5  n.m. 


RDOT:  Actual  and  noisy  range  rate  of  target  with  respect  to 

RDOTE:  interceptor,  positive  opening,  ft/sec. 


RDOTM:  Range  rate  actual,  noisy,  predicted,  and  smoothed, 

RDOTME:  n.m./sec. 

RDOTMP 

RDOTMS 


RDDOTE:  Target  noisy  acceleration,  n.m. /sec    . 


RN1 
RN2 
RN3 


Random  accelerations  of  target  for  velocity,  heading, 


and  pitch. 
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RNMP: 


RNMS: 


>  - 


Predicted  and  smoothed  range  to  target,  n0m 


RSMO: 
T: 
TA: 
TC: 

TD:* 

TM:* 
TS: 

TT:* 
VO: 

VT: 
VTDOT: 


Smoothed  range  error  from  actual,  n.m. 

Time  since  problem  start. 

Total  target  iteration  time  since  problem  start. 

Iteration  time  which  determines  time  for  plotting 

desired  limit  points. 

Time  to  target  dive . 

Target  Mach  number. 

Check  time  against  iteration  time  to  enter  main  program 

from  target  iteration. 

Time  to  target  turn. 

Intercepter  velocity ,  ft/sec. 

Target  velocity  and  acceleration,  ft/sec.  and  ft /sec    . 


VTF: 


Final  target  velocity,  a  function  of  the  target's  final 


altitude. 


*    Read  into  the  program 
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XI: 
X2: 


Points  which  mark  initial  target  location  on  the  smoothed 
plot. 


XO: 


XONM: 


XT ERR: 


Interceptor  position  along  the  north  axis,  ft.  and  n.m 


Difference  between  target's  noisy  and  smoothed  position, 


n.m. 


XTO:* 

XT: 

XTE: 

XTONM: 

XTNME: 

XTNMS: 


Target  position  along  north  axis,  initial,  actual,  and  noisy 
in  feet,  and  initial,  noisy,  and  smoothed  position,  n.m. 


XTSMO: 


Target's  difference  between  actual  and  smoothed  x  position 


Yl: 


Y2: 


Similar  to  XI  and  X2 


*    Read  into  the  program 
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YTERR: 


Similar  to  XTERR. 


YTO:* 

YT: 

YTE: 

YTONM: 

YTNME: 

YTNMS: 


Target's  position  along  the  east  axis,  similar  to  X  positions 


YTSMO:  Target's  difference  between  actual  Y  position  and  smoothed 

Y  position. 


Zl: 
Z2: 

ZETA: 


Similar  to  XI  andX2. 

Target's  deflection  angle  from  the  horizontal,  measured 
positive  downward  from  the  Intercepter,  rad. 


ZO:* 
ZONM: 


Interceptor's  altitude,  ft.  and  n.m. 


ZTERR: 


Similar  to  XTERR, 


*    Read  into  the  program 
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ZTO:* 

ZT: 

ZTE: 

ZTNM: 

ZTNME: 

ZTNMS: 

ZTSMO: 


Target's  altitude  similar  to  X  positions. 


Target's  difference  between  actual  Z  position  and  smoothed 


Z  position. 


ZTF:' 


Final  target  altitude,  ft. 


ZTPO: 


Target  pull  out  altitude,  ft 


Read  into  the  program 
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APPENDIX  B 


COMPUTER  PROGRAM  DETAILS 


A.    Purpose:    The  program  is  designed  to  generate  three  dimensional 
target  interceptor  geometries  for  use  in  analyzing  optimum  smoothing 
schemes . 


B.    Limitations:    The  interceptor  starts  at  the  coordinate  axis  origin  and 
may  have  any  northerly  heading.    The  target  may  be  placed  anywhere  in 
the  northern  semi-circle  with  any  initial  heading.    The  positive  axis 
directions  are:    X,  north,  Y,  east,  and  Z,  down,  measured  from  the 
interceptor's  altitude.    The  problem  will  stop  anytime  the  target  flies 
out  of  a  65  deg.  cone  from  the  interceptor's  heading. 


C.    Program  Capabilities:    The  program  will  compute  the  actual  range , 

range-rate,  antenna  azimuth  and  elevation  for  a  radar  frame  time  of 

two  seconds.    It  will  put  noise  on  the  above  having  the  standard 

deviations  of  2.5  n.m.  for  range,  8.72kts  for  range -rate,  CL745  deg. 

for  azimuth,  and  0.86  deg.  for  elevation.    It  also  puts  a  ha  Iff  g  acceleration 

noise  on  the  target.    It  will  compute  the  noitey  target  location,  based  on  the 
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noise.    It  will  further  transfer  the  actual  and  noisy  radar  parameters 
into  a  moving  north,  east,  down  coordinate  system  which  is  fixed  to 
move  with  the  interceptor.    It  will  determine  whether  detection  was  made 
or  not  each  frame  time,  based  on  a  detection  probability  curve,  and 
perform  a  quality  counter  operation  which  prevents  non-detected  targets 
from  being  smoothed  after  they  are  lost  for  a  reasonable  number  of  frames 

D.    Input  data  requirements:    Foiur  input  cards  are  required. 

1.  The  first  card  reads  in  16  Hollerith  characters  which  identify  the 
program  on  the  print  out  and  graphs .    The  format  is  2A8 . 

2.  The  second  card  reads  in  the  interceptor's  altitude  in  feet,  the 
target's  X  and  Y  coordinates  in  nautical  miles ,  and  altitude  in  feet, 
and  interceptor  and  target  heading  in  degrees.    The  format  is  7F6.0. 

3.  The  third  card  reads  in  the  target  and  the  interceptor's  Mach  number, 
the  target's  maneuver  g  load  capabilities  for  turn,  pull  out,  and  dive, 
(negative  number) ,  and  the  problem  time  to  target  turn  and/or  dive . 

The  format  is  7F6.2. 

4.  The  fourth  card  reads  in  a  fixed  point  number  that  allows  a  different 
starting  point  for  the  subroutine  RNDEV,  and  three  fixed  point  ones  or 
zeroes  which  control  the  graph  output.    The  first,  if  zero,  suppresses 
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the  profile  curves,  the  second,  if  zero,  suppresses  the  smoothed 
difference  plots„and  the  third,  if  zei>o,  suppresses  the  three  dimensional 
pms  smoothed  difference  plot.    The  format  is  416. 

E.  Requirements  for  reprogramming  the  smoothing  section:    The  entering 
arguments  of  the  smoothing  and  prediction  section  are  RNME ,  DIRNE , 
DIREE,  DIRDE,  and  NQC .    The  suggested  statement  numbers  to  use  are 
the  600  series  or  above  2000. 

F.  Quality  indicators:    The  primary  quality  indicator  is  the  three 
dimension  rms  smoothed  error  from  the  actual  position.    The  three 
dimensional  rms  noise  error  is  computed  for  comparison.    The  desired 
smoothing  error  on  the  individual  quantities  is  0.5  n.m.    This  could 
give  a  three  dimensional  rms  smoothed  error  of  as  much  as  0.875  n.m. 
The  output  curves  have  points  plotted  indicating  these  limits. 

G.  Output:    The  print  out  contains  the  run  identification  read  in,  the 
problem  time,  the  three  dimensional  tfms  actual  and  smoothed  error, 
whether  detection  was  made  or  not,  ND  equals  one  or  zero,  and  the 
status  of  the  quality  counter,  NQ.    There  are  nine  output  curves.    They 
plot  the  following: 
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1.  The  interceptor  and  target's  actual  and  noisy  positions  in  the 
north-east  plane. 

2.  The  interceptor  position  and  the  target's  smoothed  position  in  the 
north-east  plane. 

3.  The  interceptor  and  target's  actual  and  noisy  positions  in  the  north- 
down  plane. 

4.  The  interceptor's  position  and  the  smoothed  target's  position  in  the 
north-down  plane. 

5.  The  range  smoothed  differences  from  actual  and  the  desired  limit 
points ,  versus  time  . 

6.  The  north  direction  cosine  smoothed  differences  from  actual  and  the 
desired  limits  points  plotted  as  a  function  of  actual  range,  versus  time. 

7.  The  east  direction  cosine  smoothed  differences  similar  to  the  north 
direction  cosine  curve. 

8.  The  down  direction  cosine  smoothed  differences  similar  to  the  north 
direction  cosine  curve. 

9.  The  three  dimensional  r;ms  smoothed  error  from  actual  and  the 
desired  limits  points,  vertsus  time. 

H .    Run  time:    The  average  run  time  is  between  2  minutes  and  30  seconds 


and  2  minutes 


and 5 50!  seconds. 


»  .    i     > 

i         i    1 1    >    i     , 


»1>3         li>         •>>>  1)4  •  >    I 
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PROGRAM  THREE  D 


C 

c 
c 
c 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


THIS  PROGRAM  COMPUTES  THREE  DIMENSIONAL  I NTERCEPTER- TARGET 
GEOMkJTRItR'fOK  THE1  'PUtfPOSg'  V7\  t)i»T»?Ml!2  I  N'Q  '  V'AHOUS  PARAMETERS, 
THi  PRD(iv?^M  WiUL  fc.TOP  'ANY'T.IME  WAT  :THg'  'TA;RGET  FLIES  OUTSIDE  OF  A 
65  DEGREE  CO'n£'  FROM  '  T'HE  1'nTERC'EPTER5'  'HEaQ  TNG  , 


FOUR  INPUT  CARDS  ARE  REQUI 

1.  THE  FIRST  CARD  READS  I 
IFY  THE  PROGRAM  ON  THE  PRI 

2.  THE  SECOND  CARD  READS 
TARGETS  X  AND  Y  COORDINATE 
IN  FEET,  AND  INTERCEPTER  A 
FORMAT  IS  7F6.0. 

3.  THE  THIRD  CARD  READS  IN 
TARGET  MANEUVER  G  LOAD  CAP 
(NEGATIVE  NUMBER),  AND  THE 
THE  FORMAT  IS  7F6.2. 

4.  THE  FOURTH  CARD  READS 

A  DIFFERENT  STARTING  POINT 
FIXED  POINT  ONES  OR  ZEROS 
FIRST  IF  ZERO  SUPPRESSES  T 
SUPPRESSES  THE  SMOOTHED  DI 
SUPPRESSES  THE  3D  RMS  SMOO 


RED 

N  16  HOLLERITH  CHARACTERS  WHICH  IDENT- 

NT  OUT  AND  GRAPHS.   THE  FORMAT  IS  2A8. 

IN  INTERCEPTERS  ALTITUDE  IN  FEET, 

S  IN  NAUTICAL  MILES  AND  Z  COORDINATE 

ND  TARGET  HEADING  IN  DEGREES.   THE 

TARGET  INTERCEPTER  MACH  NUMBER, 
ABILITIES  FOR  TURN  PULL  OUT  AND  DIVE 
PROBLEM  TIME  TO  TARGET  TURN  AND  DIVE. 

IN  A  FIXED  POINT  NUMBER  THAT  ALLOWS 

FOR  THE  SUBROUTINE  RNDEV,  AND  THREE 
WHICH  CONTROL  THE  GRAPH  OUTPUT.   THE 
HE  PROFILE  CURVES,  THE  SECOND  IF  ZERO 
FFERENCE  PLOTS,  AND  THE  THIRD  IF  ZERO 
THED  ERROR  PLOT.   THE  FORMAT  IS  416, 


DIMENSION  AND  READ  STATEMENTS 


DIM 

l.ZT 

2RNM 

3YTN 

4DIR 

5,DI 

6,  DI 

7(15 

B(i5 

9D, 

DIM 

1,XT 

DO 

REA 

100  FOR 

REA 

ITT, 

11  FOR 


ENSION 

NM(151 
S(l5l) 
MS(l5l 
NE(l5l 
RNS(l5 
RNDD(1 
1),EPS 
1>,ZET 
AERRPC 
ENSION 
ERR(l5 
2  J  =  l, 
D  10  0  C 
MAT(2A 
D  11, Z 
TD,NOI 
MAT(7F 


T<151> 
),RNM(1 
,RDOTM( 
),ZTNMS 
),DIREE 
1),DIRE 
5i),DlR 
ID(l5l) 
A(l5l), 
3o),AER 

XK2), 
1)#YTER 
6 

ITU), I 
8) 

Q,XTQNM 
SE,NGRA 
6.0,/, 7 


,XONM 
51), A 
l5l>. 
(151) 
(l5l) 
S(l5i 
EDD(1 
,EPSI 
NQ(l5 
RM(3q 
Yl(2) 
R(l5l 

=  5,6) 


(151),Y0NM(151) ,Z0NM(151),XTNM(151),YTNM<151) 

TAD(151),RNME(151)»RD0TME(15D*  ATADE(151), 
XTNME(l5l  ),  YTNME(l5D,ZTNME(l5l),XTNMS(l5l), 

lT(l2),DIRN(l5i),DIRE(l5l),DIRD(l5l), 
,DIRDE(l5l),RDD0TE(l5l),RNMP(l5i),RD0TMP(l5i) 
),DlRDS(l5l),DIRND(l5l),DIRED(l5l),DIRDD(l5l) 
5l),DlRDDD(l5i),DIRNP(l5i),DlREP(i5i),DIRDP 
DE(l5l),RSM0(l5i),ANSM0(l5l),AESM0(l5l),ADSM0 
D,XTSM0(l5i),YTSM0(l5l),ZTSM0(l5i),CIRERR(l5 
),RERRP(3o),RERRM(30),TC(3o)»RERR(l5i) 

,Z1(2),X2(2),Y2(2),Z2(2),CERRP(30) 
)#ZTERR(l5i),CIRSM0(l5i),ND(l5l)   - 


,YToNM,ZTo,ZTF,PSIOD,PSITD,TM,OM,GT,GDP,GDN, 
F1,NGRAF2,NGRAF3 

F6.2,/,4I6) 


INITIAL  CONDITIONS 


300  T(1)=0. 
TC(1)=0. 
TA=0.     »   - 
TS  =  0. 

NERMAX  =  3o'1   •■ 
PI=3. 14159265 
RAU=180./PI 
DIVE=-6q./RAD 

DT=.05 

NUMPTS=l5i 

NPRINT=68 


NQC 
NQ( 
NQ( 
NCT 
F  =  l 
Q  =  i 

» 

K  =  l 
XTO 

YTO 
XTs 
YT  = 
ZT  = 
PSI 
PSI 
PSI 
AT  = 
ATD 
ATD 
BT  = 
BTD 
BTD 
RN1 
RN2 
RN3 
NUN 


=  1 

1)  =  0 
0)  =  0 

c=o 


111         >»»      ■»»> 


•  >  I 


?XT"ONM*JV7iVi         

=YTQNM*6076,1 

XTO 

YTO 

ZTO 

To=PSITD/RAO 

0=PSIOD/RAD 

T=PSIT0 

PSITO 
OT  =  0, 

DOT=0. 
0. 

OT  =  Q. 
DOT=0. 
=  0. 
=  0. 

=  0. 
IF=122q7o3125 

TARGET  AND  INTERCEPTER  VELOCITY  COMPUTATIONS 


VTDOT=Q, 

IF  C  ZTO-37000 . )63,69,69 
63  VT=TM*(lll6.-,oo398*ZTO) 

GO  TO  57 
69  VT=TM*967,6 

57  IF(ZO-37000. )73,74,74 

73  V0  =  0MM1116,-.QQ398*Z0) 

GO  TO  75 

74  V0=0M*967.6 

75  IF(ZTF~10000. >78,79,79 

78  VTF=1340.+Q.08*ZTF 
DELTAT=(VT-VTF)/56f4 

GO  TO  58 

79  VTF=VT 

DELTAT=0. 

58  CONTINUE 

C  THIS  ALLOWS  A  VARIATION  IN  THE  RANDOM  NUMBER  SEQUENCE 

IF<NOISEUo9,lQ9,l22 
122    DO    102    L=l# NOISE 

'   102  CALL  RNDEV(NUNJF,DEV) 
109  CONTINUE 

C  BEGIN  MAIN  PROGRAM  DO  LOOP 

DO  2q  I=le,4.5'l  •-  ••*    •  ";"  ,*" '  ,  •"•  •,,•  ;  *""  ,  ]     '; 

C  CHECK'  ITERAT'IOW'  T  I>iE    A'GATUST    RADAR    fRAME    T|ME 

22    IF (ABSF(TS-T A) -10. E-5) 25,25,24 

C  CHECK  MANEUVER  TIMES  AND/OR  READ  IN  MANEUVER  CONDITIONS 

62  of  10  0 


•  •  •        >  I  »     » 

t 
»  >  1  1 

>  1 

»  >         >  »  >         I 


1  I  »      1 


>  >  >      <  I  > 

1  I  I 

I 

>  >  I      III 


24    TA=TA+DT 

IF(TA-TT)39, 28/28 

39  ATDOT=0. 
GO    TO    29 

28  ATD0T=GT*32.l7/VT 

29  1F^Tfc-TD,M0i!4l,4i. 

40  ^TLiOT^O .:    J    ;,,• 
GO'  TO  42 '   ' 

41  IF(ZWTF)76,77,77 

76  BTD0T*(6DP-1.)*32.17/VT 
GO  TO  42 

77  BTD0T=(GDN-1. >*32.17/VT 
VTD0T=-56,4 

42  IF( (PSITO+PI )-AT)46,46,45 
46  ATDOT=0. 

AT=PSITo+Pl 
45  IF(8T-DIVE)34,34,35 

34  BTDOT=0. 
BT=DIVE 

35  IF  (F)  2qi,2q1,2q3 

20  3  ZTP0=(1,-C0SF(BT))*VT**2/((GDP-1. > *32 . 17 >*ZTF 
IF  (ZTP0-ZT)38,37,37 

37  IF(F)201, 201,200 

200  F=F-1. 
VTPO=VT 

201  BTDOT=(GDP-l. )*32.l7/VTP0 
IF(BT)38,43,43 

43  BTDOT=Q, 
BT  =  Q. 


38  IF(ZT-ZTF)205,2Q5,206 

205  ZT=ZTF 

206  IF<TA-<TD+DELTAT))5l,5l,2o2 

202  VTDOT=Q. 
VT=VTF 

ITERATION  OF  TARGETS  POSITION 

5i  VT=VT+VTD0T*DT+RN1*DT 

AT=AT+ATD0T*DT*RN2*DT**2/2, 

PSIT=AT 

BT=BT+BTD0T*DT*RN3*DT**2/2, 

COMPUTE  TARGET  POSITION 

XT=XT+VT*COSF(BT)*COSF(AT)*DT 
YT=YT+VT*COSF(BT)*SINF(AT)*DT 
ZT=ZT+VT*SINF(BT)*DT 

RETURN  TO  START  TO  CHECK  AGAINST  RADAR  LOOK  TIME 


25 


GO  TO  22 
TS=TS+2, 


COMPUTE, 


X0=V0*7( 


INTERCF.PTERS.  .POSITION,  ANI)  A.CTU.A'.  RAD.AR.  PARAMETERS 


59 


•I  )*O0SF'(?O2C-: 

YO=VO*T< I)*S1NF(PSI0) 

R=SQRTF( <XT-X0)**2*(YT~Y0)**2*<ZT-Z0)**2) 

EPSI=ASINF( (ZO-ZT)/R) 
IF (PSIO-PI )6q,59,59 

PSI0=PSI0'2.*PI 


4     ••  *  t 


4         i.  i     ' 


60  IF(YT-Y0)3l,32,32 

3i  ATA=-ACOSF< ( XT-XO ) / <R*COSF (EPS I ) ) ) -PS 1 0 

GO  TO  33 
32  ATA=ACOSF((XT-XO)/(R*COSF(EPSI)))-PSlO 


33 

26 
49 

48 
50 


15 
14 
27 


"CHECK"  IF'1  TARGET  ,I'3'  STTU 


•IN-'OAD 


AR 


•B£A 


M,  IF  NOT  STOP 


IF(ATA,-l,,l34,5)2'6,26;5,o"   "' 
IF( ATA+1. 1345)50, 49, 49 
IF(EPSI-i?i345)48,48,50 
IF(EPSI+1. 1345)50, 27, 27 
NUMPTS=( I-l> 
NERMAX=NUMPTS/10 
NPLOT=NUMPTS-NCTC 
IF(NUMPTS-68)l5,i4,i4 

NPRINT=NUMPTS 
GO  TO  1 

RDOT  =  -VT*COSF(PSIT-ACOSF( ( XT-XO ) / ( R*COSF ( EPS  1 ) ) )*PJ )*COSF(EPSI) 
l-VO*COSF( ATA)*COSF(EPSI ) 

CONVERT  ACTUAL  RADAR  PARAMETERS  TO  KTS,  NM/SEC,  DEG 


RNM(  I  )=R/6076.1 
ATAD( I )=ATA*RAD 
EPSID(I )=EPSI*RAD 
RDOTMU  )=RDOT/6076,1 
XONM( I)=XO/6o76,l 

YONM(  I)=Y0/6Q?6.1 
ZONM< I  )=ZO/6o76,l 
XTNMd  )=XT/6o76,l 
YTNM(  I  )sYT/6o76.1 
ZTNM< I )=ZT/6o76,l 

OBTAIN  TARGET  MANEUVER  NOISE 

CALL  RNDEV(NUNIF,DEV) 

RN1=16,*DEV 

CALL  RNDEV(NUNJF,DEV) 

RN2=16./VT*DEV 

CALL  RNDEV(NUNJF,DEV) 

RN3=l6,/VT*DEV 


OBTAIN  NOISY  RADAR  PARAMETERS  AND  CONVERT  TO  CONVENIENT  UNITS 


CALL  RN 
ATAE=AT 
CALL  RN 
RE=R+l5 

CALL  RN 
EPSIE=E 
CALL  RN 
RDOTE=R 

ATADEd 
EPSIDE( 
RNME(  I  ) 
RDOTME( 


DEV(NU 
A+,  013 
DEV(NU 
000  •  *D 
DEV(NU 
PSI+.O 
DEV(NU 
DOT+14 

)=ATAE 
I)=EPS 
=  R[=/63 
I;)=RD3 


NIF,DEV) 

*DEV 

NIF,0EV) 

EV 

NIF,DEV) 

15*DEV 

NIF,DEV) 

.7i74*DEV 

*RAD 
IE*RAD 
7.6tl-  •    • 

•TE/6q7'6w:i 


COMPUTE  NOISY  TARGET  SPACE  STABILIZED  POSITION 


IF(PSIO-l,)70,70,7l 

70  XTE=XO+RE*COSF(ATAE+PSIO)*COSF(EPSIE) 
YTE=YO*RE*SINF(ATAE+PSIO)*COSF(EPSIE) 


\    . 


71 
72 


GO  T 
XTE  = 

YTE  = 
ZTE  = 

XTNM 
YTNK 
Z''.'NK 


0  72 
XO+Rb* 
YO+RE* 
ZO-RE* 
EC  J  >=X 

EJHSZ 

COMPU 
AND  P 


COSF(ATAE-PSlO)*COSF(EPSIE) 

SlNF(ATAE-PSIO)*COSF(EPSIE) 

SINF(EPSIE) 

TE/6076.1 

7£/6o76il 

TE/:6Q7a,l; 


<  i  >  »  i  -i    •> 


>  >  >  i  >  ■>  >  > 

>     > 


TE  ACTUAL  DIRECTION  COSINES  FOR  INTERC6PT6R  WITH  ROLL 
ITCH  ZERO 


DIRN( I )=COSF(ATA)*COSF(PSlO)-SINF(ATA)*SINF(PSIO) 

DIRE(  I)=COS,F(ATA)*SINF(PSIO)*SINF(ATA>*COSF<PSIO) 
DIRD( I )=-SINF(EPSI  ) 

COMPUTE  NOISY  DIRECTION  COSINES 

DIRNEC I )=COSF(ATAE)*COSF(PSIO)-S1NF(ATAE)*SINF(PSIO) 

DIREE( I)=COSF(ATAE)*SlNF(PSlO)+SINF(ATAE)*COSF(PSlO) 
DIRDE(I)=-SINF(EPSIE) 

OBTAIN  DETECTION  PROBABILITY  AND  DETERMINE  IF  TARGET  WAS 
DETECTED 


CALL  RNDEV  <NUNIF,DEV> 

IF(RNM( I)-2o. )2l,23,23 
21  PD=2,36-.043*RNM(  I  ) 

GOTO  52 
23  IF  (RNM( I )-50, )36,47,47 

36  PD=1,95-.Q239*RNM(I) 

GO  TO  52 
47  IF  (RNM( I >-Bo. )65, 66,66 

65  PD  =  1.55-,016*RNM(  \) 
GO  TO  52 

66  IF  (RNMC I )-110, )62,68,68 
62  PD=.88-,0077*RNM(  I  ) 

GO  TO  52 

68  PD=.25-.0019*RNM(I) 

52  CONTINUE 
ADEV=ABSF(DEV) 
IF(PD-ADEV)53#54,54 

C  QUALITY  COUNTER  COMPUTATIONS 

53  ND( I )=0 

NQ( I )=NQ( I -1 )  »1 

IF  (NQ(I  ))25o,l5l,l52  ,. 

250  NQ(I)=0 
GO  TO  151 

54  ND(  I  )=1 

NQ(  I  )=NQ(  I-D+2 

IF  (10-NQ(I))25l,l53,l53 

251  NQ( I )'10 
GO  TO  153 

c  :    i    \\\  smootiukg  >m  prediction  secijion.:"; 

r>  » *>       i   mo       >i">       p  lit  \j  j      9  />  3      a  o  Q  i>  9  1  •  >   i  i       iii  i  i       i  i 

C  NO  DETECTION  CONDITIONS 

C  IF  QUALITY  COUNTER  IS  ZERO  PRINT  THE  FOLLOWING  AND  START  NEXT 

C  FRAME  TIME 


t        t    t    t  *    *    i  I 


4    .. 


151  CIRERR  < I )=999.9999 
CIKSMCK  I )=999.9999 
XTNMS( I )=XTNM( I > 
YTNMS(  I  )?YTNM( I ) 
ZTNMS( I )=ZTNM< I ) 
ND( I )=99 
NQC  =  1\ 
GO  TO  '^0 


C  IF  QLJALITY  COUNTER 

c  ones;   COMPUTE  ONLY 

152  RNME( I )=RNMP( I  ) 

RDOTME< I )=RDOTMP( I> 
DIRNE( I )=DIRNP( I  ) 
DIREE( I)=DIREP(I) 
DIRDEC I )=DIRDP( I ) 
GO  TO  6i 


IS  NON  ZERO  ERROR  VALUES  EQUAL  PREDICTED 
NEXT  PREDICTED  VALUES 

(  t 


AT  START  OR  AFTER  CONTACT  HAS 


C  DETECTION  CONDITIONS 

153  IF  (NQC)6i,6i,i54 

C  INITIAL  CONTACT  CONDITIONS 

C  BEEN  LOST  AND  REGAINED 

154  NQC=NQC-1 
RNMP( I )=RNME( I ) 
RDOTMPd  )=RDOTME(I) 
DIRNP( I )=DIRNE(  I) 
DIREP(I)=DIREE(I) 
DIKDP( I)=DIRDE(I) 
RDDOTE( I -1 ) =0 . 
DIRND( I -1 >  =0 . 
DIRNDD(I-1)=0. 
DIRED(I-1)=0. 
DIREDD(I-1)=0. 
DIRDD( I-1)=0. 
DIRDDD(I-1)=0. 

61  CONTINUE 


SMOOTHING  CONSTANT  COMPUTATIONS 

ALPHAR=.45 

ALPHAN=.43 

8ETAN=ALPHAN**2/(2.-ALPHAN) 

GAMMAN=2.*BErAN/ALPHAN-(ALPHAN+8ETAN  ) 

ALPHAE=ALPHAN 

ALPHAD=,35 

BETAD=ALPHAD**2/(2,-ALPHAD) 

GAMMAD=2.*BETAD/ALPHAD-(ALPHAD+BETAD  ) 

SMOOTHING  AND  PREDICTION  COMPUTATIONS 

R  A  N  G  E  ■>  >  '    i  »  i  o    i  •  >  o    o  o  o  i     .  ••    ■  >  >  .  ,   ,    , 

0     0  II  .  ,  J  )  I  ,  . 

*     ">     »  •>  J     0      0  0  13        •       ,  ,  ,  ,  ,       , 

RNMS(  I  )=RNMP(  ;**ALPHAfUCflNM6<  I  J.tRNMP(  I  )  ) 

RDDOTE( I )=RDDOTE( I *1 ) +1 . / ( 2 . *Q ) * ( RDOTME ( I  )-RDOTMP<  I  )  ) 

RNMP( I+1)=RNMS( I >*2.*RD0TME< I > +2 , *RDDOTE (  I  ) 
RDOTMP(  I+D=RDOTME(  I  ) +2  .  *RDDOTE  (J  ) 


NORTH  DIRECTION  COSINE 


•  •  • 

•  t 


♦  •  • 

•  * 


D I RNS ( I )=DIRNP( I )+ALPHAN*(DIRNE( I >-DIRNP< I ) ) 

DIRND(  I  )  =DIRND( I -1 ) +2 . *D I RNDD ( I -1 ) +BETAN/2 . * < DI RNE d )-DlRNPd>> 

DIRNDD(I)=DIRNDD( I -1 ) +2 . *GAMMAN/4  .  * ( D 1 RNE ( I )-DIRNP( 1)) 

DIRNP(  I*D=DIRNS(  I  >*2.*DIRND<  I  ) +2  .  *DI  RNDD  (  I  ) 


EAST  DIRECTION  CC3IN2 


■  .  i 


DIRES( I  )=DIREP( I  )*ALPHAN*(DIREE( I )-DIREP( I ) ) 

DIRED(I)=DIRED<I-1)+2.*DIREDD(1-1)*BETAN/2,MDIREE<I)-DIREP(I>> 
DIREDD( I )=DIREDD( I -1 > +2 . *GAMMAN/4 . * < D I  REE < I )-DIREP(I >> 
DIREP(  I*i) sD IRES ( I)+2.*DIRED(I)+2.*DIREDD(I ) 

DOWN  DIRECTION  COSINE 

DIRDSC I >=DIRDP( I > +ALPHAD* ( D I RDE < I )-DIRDP( I )) 

DIRDDC  I  )=DIRDD( I-1)+2.*DIRDDD< I -1 ) +8ETAD/2 . * ( D I RDE < I )-DIRDP(I>) 
DIRDDDd )=DIRDDD( I -1 ) +2 , *GAMMAD/4 . * ( DI RDE ( I )-DJRDP< I ) ) 
DIRDP(I+l)=DIRDS(I >+2.*DIRDD(I ) +2 . *D I RDDD d  ) 

COMPUTE  SMOOTHED  TARGET  POSITION 

ZETA( I )=ACOSF(DIRDS( I ) )-PI/2. 

XTNMS( I)=XONM< J)*RNMS< I ) *COSF ( ZET A ( I ) ) *D I RNS ( I ) 
YTNMS( I )=YONM( i>+RNMS< I )*COSF(ZETA< I >  >*OIRES( 1 ) 
ZTNMSC I )=ZONM(I)*DIRDS(I )*RNMS(I ) 

COMPUTE  SMOOTHED  AND  ERROR  DIFFERENCES  FROM  ACTUAL 


C 
C 


XTSMO( I )=XTNM( I )-XTNMS( I ) 
YTSMO( I )=YTNM( I >-YTNMS( I ) 
ZTSMOd  )=ZTNM(I  )-ZTNMS(I  ) 
XTERR( I  )=XTNM( I  )-XTNME( I  ) 
YTERRd  )=YTNMd)-YTNME(I  ) 
ZTERR( I)=ZTNM( I )-ZTNME<  I  ) 
RSMOCI )=RNM(I )-RNMS(I  ) 

COMPUTE  THE  ACTUAL  AND  SMOOTHED  3D  RMS  ERROR  VALUES 

CIRERR( I )=SQRTF(XTERR(I )**2  +  YTERR( I >**2*ZTERRd  >**2) 
CIRSMO( I )=SQRTF(XTSMO( I  )**2*YTSM0( I ) **2+ZTSMQ < D**2) 

CHECK  TO  SEE  IF  THE  ACOSF  ARGUMENT  HAS  BEEN  EXCEEDED* 

IF  IT  HAS  SET  THE  COSINE  EQUAL  TO  ONE  AND  PRINT  OUT  THE  VALUE 


IF  <DIRNSCI>-1.)211,211,210  ,, 

210  PRINT  212, (T( I ),DIRNS< I ) ) 

212  F0RMAT(6X2HT=  F8  .  4 , 2X6HD I RNS  =  F8 . 4 ) 

DIRNS( I ) =1 . 0 

211  IF(DIRES(I)-i,)22i, 22l,22o 

220  PRINT222,<T< I),DIRES( I)) 

222  F0RMAT(6X2HT=F8.4,2X6HDIRES=F8.4) 
DIRES(I)=l.o 

221  IF(DIRDSd)-l, 0)?3i, R33.  2?0      .   . 

230  PRINT232",  <!t  (4-5,DIRD£H  S)  '•<•   •         ;     ■' .  '  •  • 
232  F0RMAT(6X2H7*F8;4,2*XuHC<RC3  =  F3.4) 

DIRDS( I )=1.0 

231  CONTINUE 

CONVERT  THE  DIRECTION  COSINES  TO  DEGREES.  FOR  GRAPH  OUTPUT 


•  ft  « 

ft   • 


•  *  • 

•  t 


ANSMO(I)=ACOSF(DIRN(I))-ACOSF(DIRNS(I>) 
AESMO(I)=ACOSF(DIRE(I))-AC0SF(DIRES(I>) 
ADSMO< I)=ACOSF(DIRD( I ) ) -ACOSF < D I RDS ( I ) ) 

C  COMPUTE  THE  DESIRED  LIMIT  POINTS  FOR  THE  GRAPH  OUTPUTS 

IFaT(|(K/-T(  I  )  U20,l-'0>121 

120  TccK+n^rcvK'^rG'.  •  •   >j*  »•  *        .:.;,..  :.. ,• 

AERRP(K)=ATANF( ,5/RNM< I ) ) 
AERRM(K)s-AERRP(K) 
RERRP(K)5,5 
RERRM<K>=-.5 
CEKRP(K)=.875 
K  =  K  +  1 
121  CONTINUE 

C  IF  DETECTION  WAS  MISSED  INCREASE  Q  IN  THE  RDOT  SMOOTHING 

C  COMPUTATIONS 

IF  (ND( l)>55,55,56 

55  Q=Q+1. 
GO  TO  20 

56  Q=l, 

2q  T(I+l>=T(I)+2. 


END  OF  MAIN  PROGRAM  DO  LOOP 


1  CONTINUE 


OUTPUT  CONTROL 


C  COMPUTE  MARKERS  WHICH  INDICATE  THE  TARGETS  INITIAL  POSITION 

C  ON  THE  SMOOTHED  PROFILE  CURVES   - 

Xl(l)=XT0NM+.5 
Xl(2)=XT0NM-.5 
Y1(1>=YTQNM+,5 
Y1(2)=YT0NM-,5 
Zl(l)=ZTO/6o76,i*.o5 

Zl(2)=ZTo/6o76,i-.05 

X2(1)=XT0NM+.5 

X2(2)=XToNM-.5 

Y2(l)=YT0NM-.5 

Y2(2)=YToNM*.5 

Z2(l)=ZTO/6o76,i-,o5 

Z2(2)=ZTo/6o76,1+.05 

C  GRAPH  STATEMENTS 


IF  <NGRAF1)27q,270,27i 

271  CONTINUE 

DO  101  I=l»4 
101  IT(I)=8H 

IT(l)=8HGENTZ  NE 

IT(2)=8H  PLANE  A 

IT(3)=8HCTUAvLll*  >  ;  ;    ,   .. 

IT(4)=8HNCIS£.  -  -  ,,o   :.  J  . 

1=1 

CALL  DRAW(NUMPTS,YQNM.XONM,1,0#4H 
1LAST) 
CALL  DRAW(NUMPTS,YTNM,XTNM,2,0,4H 
CALL  DRAW(NUMPTS,YTNME,XTNME,3,q,4h 


, IT, l5.,l5f. 0,3,2,2,6, 6,1, 


- 


X03 


10* 


110 


270 


282 


105 


106 


116 


DO  10 
IT(  I  ) 

IT(1) 
I  T  (  2  ) 
I  T  (3  ) 

IT  (4.) 

1=1  ; 

CAl  l.k 
1LAST) 
CALL 
CALL 
CALL 

DO  10 
I  T  C  I  ) 

IT(i> 
IT(2) 
IT(3) 
ITC4) 

1=1 
CALL 
1LAST) 
CALL 
CALL 
DO  11 
IT(  I) 
I  T(  1) 
IT(2) 

I  T  (3) 
IT(4) 

1=1 
CALL 

1LAST) 
CALL 
CALL 

CALL 
CONTI 
IF  (N 
CONTI 
DO  10 
I  T  <  I  > 
IT(1) 
IT(2) 
IT  ( 3 ) 
I  T  (4  ) 

1=1 
CALL 

CALL 
CALL 

DO  10 
IT(  I) 

IT(1) 
I  T  (2) 

IT  C  3) 

IT(4) 

1=1  - 
CALL 

CALL 
CALL 

DO  11 
ITU) 

IT(1) 
IT(2) 


3  1=1,4 

=  8H 

=8HGENTZ  NE 

=8H  PLANE  S 

=8HM00THED 

=  8HPQS,IT„I0N 

r.RA  W  i  MUMPTS ,  >  ONM  '/xON*  ;  1 ,  C  *  4H     ,  !  T §£$  , ,  l5  ,  ,  0  , 3 ,  2  , 2 ,  6 ,  6  r  1 , 

\ 

DRAW(NUMPTS,YTNMS,XTNMS,2,Q,4H     ) 
DRAW  <2,Y1,X1,2,0,4H     ) 
DRAW  (2,Y2,X2,3,q,4h     ) 

4  1=1,4 

=  8H 

=8HGENTZ  ND 
=8H  PLANE  A 
=8HCTUAL  t 
=8HN01SE 


DRAW(NUMPTS,XONM,ZONM,1,0,4H 

DRAW(NUMPTS,XTNM,ZTNM,2,0,4H 
DRAW(NUMPTS,XTNME,ZTNME,3,0#4H 
0  1=1,4 

=  8H 

=8HGENTZ  ND 
=8H  PLANE  S 

=8HM00THED 
=8HP0SITI0N 

DRAW(NUMPTS,XONM,ZONM,1,0#4H 


, IT, 15, ,1.5, 0,0*2,2,6,6/1  , 


IT, 15., 1,5, 0,0*2,2,6, 6,1  , 


DRA 
DRA 

DRA 
NUE 
GRA 
NUE 
5  I 
=  8H 
=  8H 
=  8H 
=  8H 
=  8H 


W(NUMPTS,XTNMS,ZTNMS,2,0,4H 
w(2,Xl,Zl,2,0,4H     ) 

W  (2,X2,Z2,3,0,4H     ) 
F2)281,28l,282 

=  1,4 

GENTZ  RA 
NGE  SMOO 
THING  RE 
SULTS 


DRAW(NUMPTS,T,RSM0,1,Q,4H 

DRAW(NERMAX,TC,RERRP,2,2,4H 

DRAW(NERMAX,TC,RERRM,3,2,4H 

6  1=1,4 

=  8H 

=8HGENTZ  NO 
=8HRTH  SMOO 
=8HTHING  RE 
=8HSULTS 

DRAW.  NUMPTS>T;ANSMO,  1,0  •,'*♦* 

ORAWiNERKVWVTCvflBRRPifcifcvAH 
DRAW(NERMAX,TC,AERRM,3,2,4H 

6    1=1,4 
=  8H 

=8HGENTZ  EA 
=8HST  SMOOT 


, IT, 5q,,1., 3,0,2,2,6, 6,1, LAST) 
) 


,11, -*o 
)' 

) 


,•.01*3, 0,2, 2, 6, 6,1, LAST) 


108 


IT 
IT 
1  = 
CA 
CA 
Ca 
DO 
IT 
IT 
IT 
IT 
IT 
1  = 
CA 

CA 
CA 

CO 
IF 

CO 
DO 
IT 
IT 
IT 

IT 
IT 

1  = 

CA 

1LA 

CA 

283  CO 


281 


284 


285 


(3)  = 

(4)  = 

1 
LL  D 

Li;  :o 

LL  D 
'  108 
(  I)  = 

(1)  = 

(2)  = 

(3)  = 

(4)  = 

1 

LL  D 
LL  D 
LL  D 

NTIN 
(NG 

NTIN 
285 

(I)  = 

(1)  = 

(2)  = 

(3)  = 

(4)  = 

1 

LL  D 
ST) 
LL  D 
NTIN 


8HHING 
8HULTS 


RES 


RAW(NUMPTS»T,AESMO,1,0#4H 

Rmwumcpmax   re,  VHWii2;a#4fi 

&W(NEPMAx:TC,A'fcRRM,3,2,*H 

1=1*4 
8H 

8HISENTZ  DO 
8HWN  SMOOT 
SHYING  RES 
8HULTS 


RA 
RA 
RA 

UE 
RA 

UE 
I 

8H 
8H 
8H 
8H 
8H 


W(NUMPTS,T,ADSMO,1,0,4H 

W(NERMAX,TC,AERRP,2,2,4H 

W<NERMAX,TC,AERRM,3,2,4H 

F3)283,283,284 

=  1,4 

GENTZ  3D 
RMS  SMO 
OTHING  R 
ESULTS 


RAW(NUMPTS,T,CIRSM0,1,Q,4H 

RAW(NERMAX,TC,CERRP,3,2,4H 
UE 


,  I  T,5o,,. 01* 3,0 #2,2,6,6,1* LAST) 
) 
)  '     : 


>  IT, 50., .01,3, 0,2, 2,6, 6,1, LAST) 
) 
) 


,IT,5o., .5,0,0*2,2,6,6,1, 


95 


64 


107 


PRJNT  STATEMENTS 

PRINT  95 

FORMAT (1H1, 27X1 7HSM00TH I NG  RESULTS) 

PRINT  64, ( IT( I), 1=5,6) 

F0RMAT(/,32X,?A8) 

PRINT  107 

FORMAT(   /,6X66HT=TIME,  R  NM=ACTUAL 
1  ACTUAL  ERROR,      , / , 6X64HC I RSMQ=3 
2,  ND  =  1  DETECTION,  ND  =  0  N0,/,6X  68HDETECT I  ON,  NQ  =  QUALITY  COUNTER,  W 
3HEN  NO=0  TARGET  CONSIDERED  LOST  AND  ,/,*X24HALL  NINES  ARE  PRINTED. 


RANGE,  CIRERR=3DIMENSI0NAL  RMS 
DIMENSIONAL  RMS  SMOOTHED  ERROR 


4      ) 

PRINT  96 

96  FORMAT( 

12X2HNQ) 
PRINT  97, 

97  FORMAT  ( 
2  CONTINUE 

400  STOP 
END 
END 


/#l5XlHT,4x4HR  NM, 7x6hC I RERR, 4x6HC I RSMO, 2X2HND, 

(T( I ),RNM(I),CIRERR(I)#CIRSMO(I),ND(D#NQ(I)#I=l#NUMPTS) 
7X,1F10.0#3F10.4,2I4) 


* : 


•  •  •  • » » • 


«  1  «  il      >     )  K 

.1 

r)        J       <J  1         J  8 


i     I     ) 


I  >>.!.> 


M 


Appendix  C 


1      i 
i    i  >  >  >    \  i  i 


Computer  Program  of  Program  ALPHA  CHECK 

Only  the  smoothing  and  prediction  and  output  sections  are  shown 
since  the  preceeding  portions  are  the  same  as  in  THREE  D. 


I'    I 


V.7I  cf  120  OC 

t  *  I     .  Q 


*■       * 


w 


<      * 


47  IF  (RNM( I  )-80. )65,66,66 

65  PD=1.55-.016*RNM(  I  ) 
GO  TO  52 

66  IF  (RNM(  I  >-110 , )62,66,68 

62  PDs.88-.0077*RNM(  I  ) 
GO  -TO  52  , 

68  PD=.25.-it&0l9*RNH(  l) 

52  CONTINUE 
ADEV=ABSF(DEV) 

ADEV=0. 
IF(PD-ADEV)53,54,54 

53  NQ(I )=0 
GO  TO  20 

54  NQ(I)=1 

20  T(  I+1)=T(  I)*2. 
1  CONTINUE 


.i  >  1 1 


SMOOTHING   COMPUTATIONS 


FPTS=NUMPTS 
NALPH=20 


302 
304 


303 
308 


3o6 


55 

56 

305 

301 


114 


604 


115 


RANGE  SMOOTHING  CONSTANT  ANALYSIS 

ALPHAR=,o5 

DO  301  J=1#NALPH 

P  =  l. 

AR( J)=ALPHAR 

ERAR( J)=o  . 

DO  3q5  i=1,NUMPTS 

IF  (NQ( I ) )3o2,302,3o3 

IF(P)304, 304,305 

RNMEd  )=RNMP<  I) 

RDOTME( I )=RDOTMP( I ) 

GO  TO  3o6 

IF(P)3o6,3o6,3q8 

P=P-1. 
CTC=T(I)/2.*1. 

RNMP(  I  )=RNME(  I  ) 
RDOTMP( I )=RDOTME(  I  ) 

RDDOTE(I-1)=0. 

RDDOTECI )=RDDOTE(I-l)+l./(2.*Q)*(RDOTM£( I )-RDOTMP( I ) )  - 

RNMS( I )=RNMP( I )+ALPHAR*(RNME( I )-RNMP(  I  ) > 

RNMP( I*l)sRNMS< I )*2.*RDOTME< I )+2,*RDD0TE< I ) 

RDOTMP( I+l)=RDOTME< I )*2.*RDD0TE(  I) 

RSMO( I)=(RNM( I )-RNMS( I  )  )**2 

ERAR( J)=ERAR( J)*RSMO( I )/(FPTS-CTC) 

IF(NQ(I))  55,55,56 

Q=Q*1. 

GO  TO  3o5 

Q  =  l. 

CONTINUE 

ALPHAR=ALPHAR+,05 
1=1 

PRINT  11 4  .  . 

FCRMrATU^i^l7y,o6K0PT,IMUM,  SMOOTHING  CONSTANT  ANALYS  I  S, //,  6X62H 
lTY  {.NDJCATC'R^T-HE  MINIMUM  AVERAGE  MEAN  SQUARE  ERROR 
PRINT  604, < < IT< I ), I=5,6),AN0ISE) 
FORMAT(//,2  0X,2A8,2X7HANOISE=F4.2) 
PRINT  115 
F0RMAT(//,29Xl6HRANGE,  ALPHAR      )  4 

PRINT  116 , <AR(J),ERAR< J), J=1,NALPH) 


QUAL 

) 


I  .  .         II.      Ill      III. 


/-irr uivuls*.  k~> 


116  FORMAT(//, ( (23X7HALPHAR=  FA  .  2 , 2X5HERAR=   F6.4)/)) 

C  NORTH  DIRECTION  COSINE  SMOOTHING  CONSTANT  ANALYSIS 

ALPHAN=.05 

DO    3  2*,.  ItsjLvNALPH.  " 

anjkj-aiphm.m  ••-     

ERAN(K)=0. 

BETAN=AlPHAN*«2/(2.-ALPHAN> 

GAMMAN=2.*BETAN/ALPHAN-(ALPHAN+BETAN  ) 

DO  325  I=i,NUMPTS. 

IF(NQ( I ) )322, 322,323 

IF(P>324,324,325 

DIRNEd  )=DIRNP(I  ) 

GO  TO  326 

IF(P>326,326,328 

P=P-1. 

DIRNP( I )=DIRNE(I) 

DIRND(I-D=0. 

DIRNDD( I-1)=0. 

DIRND(I )=DIRND(I-1)+2.*DIRNDD(I-1)+BETAN/2.*(DIRNE(I)-DIRNP(I)> 

DIRNDD( I )=DIRNDD( I -1 ) +2 . *G AMMAN/4 . * ( DI RNE ( I ) -DIRNP ( I ) ) 

DIRNS( I )=DIRNP( I )♦ ALPHA N*(DIRNE( I )-DlRNP( I )) 

DIRNP(  I+l)=DIRNS( I)+2.*DIRND( I )*2.*DIRNDD( I > 

DNSMOC I)=(DIRN(I)-DIRNS(I))**2 

ERAN(K)=ERAN(K)*DNSMO(I )/(FPTS-CTC) 

CONTINUE 

ALPHAN=ALPHAN*,05 

1=1 


322 
324 

323 

328 


326 


325 

321 

117 
118 


) 


332 
334 

333 
338 


336 


335 
331 


PRINT  117 

FORMAT(//,20X3lHNORTH  DIRECTION  COSINE,  ALPHAN 

PRINT  118* (AN<K),ERAN(K),K=1,NALPH) 

FORMAT(//, <<23X7HALPHAN=  F4 . 2, 2X5HERAN=  F1Q.8)/)) 

EAST  DIRECTION  COSINE  SMOOTHING  CONSTANT  ANALYSIS 

ALPHAE=,05 

DO  331  L=l,NAlPH 

P  =  l. 

AE(L)=ALPHAE 

BETAE=ALPHAE**2/(2,-ALPHAE) 

GAMMAE=2.*BETAE/ALPHAE-(ALPHAE+BETAE) 

ERAE(L)=0 , 

DO  335  I=1,NUMPTS 

IF  (NQ(I ) )332,332,333 

IF(P>334,334,335 

DIREP(I)=DIREE(I) 

GO  TO  336 

IF(P)336,336,338 

DIREP(I)=DIREE(I) 

DIRED(  1  -1 )  =  0  . 

DIREDD(I-1)=0. 

DIREDr(  I  /-JDI,aED(  J-1;+2.-D;^cCD(  I  -li+£ETAE/2.  >€DIREE(  I  )-DlREP(I)>- 

DIRED'0(  IV=DlREDT'(l-l>;+?.*G4WE/4  * ( Di PEE ( O -DIREP < J  )) 

DlRfcS(l)l=DlREP(l  )*aLPHAE»(DIREE(  I  /-OiREP<  1/  ) 

DIREP( 1*1 ) =DIRES( I >+2.*DIRED( I ) +2  ,  *DI  REDDU  ) 

DESMOC J )=(DIRE( I )-DIRES( I ) )**2 

ERAE(L)==ERAE(L)+DESMO(I)/(FPTS-CTC) 

CONTINUE 

ALPHAE=ALPHAE*,05 


•  •  • 
* 
•  t 

I 

4    •    • 


119 


161 


342 
344 

343 
348 


346 


345 
341 


121 
123 
271 
110 


AFPENDIX  C 
PRINT  H9 

FORMAT(1H1,20X31H  EAST  DIRECTION  COSINE,  ALPHAE 
PRINT  161,  (AE(L),ERAE(L),L  =  1,NAI_PH) 
F0RMAT(//,((23X7HALPHAE=  F4  .  2, 2X5HERAE  =  F10.8)/)) 


111 


112 


I  DOWi^.lD.IREC.TlON^ClDS'^E'eMODJH^NG  (?O^STANT  ^ANALYSIS 

ALPHAD=.05 

00  341  M=1,NALPH 
P  =  l. 
AD(M)=ALPHAD 

BElAD=ALPHAD**2/(2t.ALPHAD) 

GAMMAD=2.*BETAD/ALPHAD-(ALPHAD*BETAD  ) 
ERAD (M)=o, 

DO  345  l=i,NUMPTS 
IF(NQ( I ) )342,342,343 
IF(P)344,344,345 
DIRDEC I )=DIRDP(1) 
GO  TO  346 

IF(P)346,346,348 

P=P-1. 

DIRDP( I )=DIRDE(  I) 

DIRDD( I-1)=0. 

DIRDDD(I-1)=0, 

DIRDD<I>=DIRDD(I-l>+2.*DlRDDD(I-l)*BETAD/2fMDIRDE(I)-DlRDP<I>> 
DIRDDD(I)=DIRDDD(  I-1)*2.*GAMMA0/4.*(DIRDE(I  )-DIRDP(D) 
DIRDS(I)=DIRDP( I >+ALPHAD*(DlRDE( I ) -D I RDP( I ) ) 

DIRDP(I*1)=DIRDS(I)*2.*DIRDD(I)*2.*DIRDDD(I) 

DDSM0(I)=(DIRD(I)-DIRDS(I))**2 

ERAD(M)=ERAD(M)+DDSMO(I )/<FPTS-CTC) 

CONTINUE 

ALPHAD=ALPHAD*,05 

1=1 

PRINT  121 

FORMAT(//,20X3lHD0WN  DIRECTION  COSINE,  AL.PHAD   ) 
PRINT  123,(AD(M),ERAD(M),M=1,NALPH) 
F0RMAT(//,((23X7HALPHADs  F4 . 2, 2X5HERAD=  F1Q.8)/)) 
IF  <NGRAPH)27o,27o,27l  * 

CONTINUE 

DO  110  1=1,4 
IT( I )=8H 

IT(1)=8HGENTZ  AL 
IT(2)=8HPHAR  VS 
IT(3)=8HAVG  MEAN 
IT(4)=8HSQ  ERROR 
1=1 

CALL  DRAW<NALPH,AR,ERAR,Q,Q,4H     , I T, . 2, 0 . , Q. 0 i 2, 2, 5, 6, 1 , LAST ) 

DO  111  1=1,4 

IT(I)=8H 

IT(1)=8HGENTZ  AL 
IT(2)=8HPHAN  VS 
IT(3)=8HAVG  MEAN 
IT(4)=8HSQ  ERROR 

1=1      

CALL  DR.WrNALPH.AN.SRA'M^O,^     ,  I  T  ,  .  2,  C  . ,«  ,>*,  2i  2,  5,  6, 1 ,  LAST  ) 

I T( I >=8H 

IT(1)=8HGENTZ  AL 
IT(2)=8HPHAE  VS 
IT(3)=8HAVG  MEAN 
IT(4)s8HSQ  ERROR 


«    f  9    *    f 


t    *    t  ft 


CONFID 


1=1 


APPENDIX  C 


it  i»i 

i       t 

i  * 


•    *    »       »    >  *    «    *    >       » 


lit       • 


CALL    DRAW(NALPH,AE,ERAE,Q#Q'4H 

DO    113    1=1,4 
113    IK  I  )-.8H 

iTiii^npfiNTziAL  :,,: 

ITC2)i,8HPHAO    VS  '    * 
IT(3)=8HAVG    MEAN 
IT(4)=8HSQ    ERROR 

1=1 

CALL    DRAW(NALPH,AD,ERAQ,o,Q,4H 

270    CONTINUE 

AN0ISE=AN0ISE*,5 

2    CONTINUE 
END 
END 


,IT,.2,Q.,o>0»2,2,5,6,i,LAST) 


>  i   1)9 


,  IT,  .2,0. ,0,0*2, 2, 5, 6,1, LAST) 


j  o 
a  ft 


*   <A 


#    #    *  +    +    + 


Jit    »  »   »  9       9 

(I  0  9  ii 


»    .  . i    '       i< 


•        i  ,i    «    «        »    »    i 


»•  ",'    '     '    "  'Append IX' C 


iii       •   i   i   i       iii 


Results  of  Runs  2A  and  3A 


JO)       01 


7£.  of,  Ldn 


I         ©     >  O  J  X 


1  .'  y    ■'       J  »    I       J    t 


•  -J  *      >  x> 


C  (  I  ... 


•  •  t 
•  1 


±&£>1jU     V  J. 


SMOOTHING  RESULTS 


RUN  2A 


T=TIME,    R    NM  =  ACTUAL    RANGE,    C I RERRs3D I  MENS  I ONAL    RMS    ACTUAL    ERROR, 
CIRSM0*3  .DIMENSIONAL    RWS    SN00TH€9.  -ESR-ORv  -ND-l  -DSTECT  I  ON,    ND»0    NO 

jN&tfQUAl;  I  TY-j COUNTER;  •  ,WH£N ;  N-Q  =  Q;  TA-flGB-T    CONSIDERED    LOST    AND 

ARE  >S'If1T£'D:.  •  ••■  »' : •", 


DETECTjOH, 
ALL  NINE'S 


2, 

1     14, 

4, 

14, 

6, 

15, 

8, 

15. 

10. 

16, 

12, 

16, 

14. 

17. 

16, 

17. 

18, 

17, 

20, 

18, 

22, 

18. 

24, 

19. 

26, 

19. 

28, 

20. 

30. 

20. 

32, 

21. 

34, 

21. 

36, 

22, 

38. 

22, 

«o. 

23, 

42, 

23, 

44, 

24, 

46, 

24, 

48, 

25, 

50. 

25. 

52, 

26. 

54, 

26, 

56, 

27. 

58, 

27, 

60. 

28, 

62, 

29, 

64, 

29. 

66. 

30, 

68, 

30. 

70, 

31, 

72. 

31, 

74, 

32, 

76. 

32, 

78, 

33. 

80, 

33. 

82, 

34, 

84, 

34. 

86, 

,:  '53; 

88; 

3' 5'. 

90, 

'  3S°. 

92. 

36. 

94, 

36. 

96, 

37, 

98, 

37, 

100, 

38. 

R  NM 
1508 
5361 
9288 
3280 
7357 

1530 
5815 

0281 
4867 
9446 
4l55 
8975 
3867 
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